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Abstract 

Experimental  and  modelling  studies  have  been  conducted  on  catalytic  autothermal  reforming  (ATR)  of  methane  for  hydrogen  production  over 
a  sulfide  nickel  catalyst  on  a  gamma  alumina  support.  The  experiments  are  performed  with  different  feedstock  under  thermally  neutral  conditions. 
The  results  show  that  the  performance  of  the  reformer  is  dependent  on  the  molar  air-to-fuel  ratio  (A/F),  the  molar  water-to-fuel  ratio  (W/F)  and 
the  flowrate  of  the  feedstock  mixture.  The  optimum  conditions  for  high  methane  conversion  and  high  hydrogen  yield  are  A/F  =  3-3.5,  W/F  =  2-2.5 
and  a  fuel  flowrate  below  120-250 1  h_1 .  Under  these  conditions,  a  methane  conversion  of  95-99%  and  a  hydrogen  yield  of  39-41%  on  a  dry  basis 
can  be  achieved  and  1  mole  of  methane  can  produce  1.8  moles  of  hydrogen  at  an  equilibrium  reactor  temperature  of  not  exceeding  850  °C. 

A  two-dimensional  reactor  model  is  developed  to  simulate  the  conversion  behaviour  of  the  reactor  for  further  study  of  the  reforming  process. 
The  model  includes  all  aspects  of  the  major  chemical  kinetics  and  the  heat  and  mass  transfer  phenomena  in  the  reactor.  The  predicted  results  are 
successfully  validated  with  experimental  data. 

©  2005  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Fuel  cell  technology  has  been  of  special  interest  in  recent 
years  due  to  its  higher  efficiency  in  energy  conversion  and  no/less 
harmful  emissions  compared  with  other  energy  conversion  sys¬ 
tems.  All  fuel  cells  are  either  fed  by  pure  hydrogen  from  a  storage 
vessel  or  by  hydrogen-rich  gas  produced  from  widely  available 
fuels  via  a  reforming  process.  A  direct  supply  of  pure  hydrogen 
is  desirable  for  fuel  cell.  Nevertheless,  hydrogen  storage  still 
remains  a  challenge  as  the  fuel  has  a  very  low  energy  density 
under  normal  ambient  conditions,  and  this  makes  storage  dif¬ 
ficult  for  mobile  applications.  Therefore,  a  reliable  method  to 
ensure  a  steady  supply  of  hydrogen  for  mobile  fuel  cells  is  to 
use  reforming  techniques,  which  extract  hydrogen  from  hydro¬ 
carbon  fuels  such  as  methane. 

Three  major  thermo-chemical  reforming  techniques  are  used 
to  produce  hydrogen  from  methane  and  other  hydrocarbon  fuels 
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[1,2],  i.e.,  steam  reforming  (SR),  partial  oxidation  (PO)  and 
autothermal  reforming  (ATR).  The  steam  reforming  of  methane 
has  been  studied  extensively  [3-6].  It  is  probably  the  most 
common  and  traditional  method  for  producing  hydrogen  on  an 
industrial  scale.  Though  this  process  can  yield  high  a  concen¬ 
tration  of  hydrogen  (up  to  70%  on  a  dry  basis),  it  is  strongly 
endothermic  and,  hence,  requires  a  substantial  supply  of  exter¬ 
nal  heat.  Therefore,  a  reforming  system  with  a  heat- exchanger 
becomes  very  bulky  and  heavy,  and  it  has  high  thermal  iner¬ 
tia  for  frequent  start-up  and  shutdown  operation.  As  a  result, 
it  is  not  so  suitable  for  a  mobile  fuel  cells.  Partial  oxida¬ 
tion  [7-10]  does  not  have  the  disadvantage  of  being  endother¬ 
mic,  but  it  produces  a  high  carbon  monoxide  concentration 
[7]  that  is  undesirable  for  polymer  electrolyte  membrane  fuel 
cells. 

Autothermal  reforming  [11-16]  combines  the  thermal  effects 
of  the  PO  and  SR  reactions  by  feeding  the  fuel,  water  and  air 
together  into  the  reactor.  The  thermal  energy  generated  from  PO 
is  absorbed  by  SR  and  hence  the  overall  temperature  is  lower. 
This  is  favourable  for  the  water-gas  shift  reaction  which  con¬ 
sumes  carbon  monoxide  and  produces  more  hydrogen  [7,11]. 
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Nomenclature 

b 

thickness  of  the  reactor  wall  (m) 

Ci 

mole  fraction  of  gas  (/) 

cpg 

heat  capacity  of  gas  (Jkg-1  K-1) 

cpb 

heat  capacity  of  catalyst  bed  (Jkg-1  K-1) 

Di 

gas  diffusivity  of  species  i  (m2  s-1) 

^dp  i 

dispersion  coefficient  (m2  s-1) 

dip 

catalyst  pellet  diameter  (m) 

Ej 

activation  energy  of  reaction  j  (kJkmol-1) 

A  Hi 

adsorption  enthalpy  of  species  i  (kJkmol-1) 

A  Hj 

heat  of  reaction  j  (kJkmol-1) 

h 

heat  transfer  coefficient  (W  m-2  K- 1 ) 

h,  h0 

heat  transfer  coefficient  on  the  inside  and  outside 
of  reactor  wall  (W  m-2  K-1) 

K 

heat  conduction  coefficient  of  catalyst  bed 
(W  m-1  K_1) 

Kej 

equilibrium  constant  of  reaction  j  (j  =  2,  3,  4) 

Kcu4 

adsorption  constant 

Koi,  K0j 

constant 

*ch4 

adsorption  constant 

kg 

heat  conduction  coefficient  of  gas  (W  m-1  K-1) 

ki 

rate  constant  of  reaction  j 

k0j 

constant 

PC  h4 

partial  pressure  CH4  (bar) 

R 

universal  gas  constant  (kJkmol-1  K-1) 

Rj 

rate  of  reaction  j  (kmolkgcat-1  h-1) 

r 

radial  coordinate  (m) 

n 

conversion  rates  of  species  i  (kmolkgcat-1  h-1) 

Sh 

heat  transfer  area  per  volume  of  catalyst  bed 
(m2  m-3) 

T,Tg 

temperatures  of  catalyst  and  bulk  gas  (K) 

Ta 

surrounding  temperature  (K) 

Ta 

outer  reactor  wall  and  surrounding  temperatures 
(K) 

u 

superficial  gas  velocity  (ms-1) 

Greek  symbols 

a 

overall  heat  transfer  coefficient  through  the  reac¬ 
tor  wall  (W  m-2  K-1) 

£ 

void  fraction  of  catalyst  bed 

X 

heat  conduction  coefficient  through  the  wall 

Pc  at 

catalyst  density  (kgm-3) 

Pb 

catalyst  bed  density  (kgm-3) 

Pg 

gas  density  (kgm-3) 

Subscripts 

a 

ambient  air 

cat 

catalyst 

g 

gas 

i 

gas  species 

j 

reaction  index  (1-4) 

r,  z 

axes  cylindrical  coordinate 

s 

solid  phase 

Hence,  the  autothermal  reactor  is  more  compact  and  practical 
for  use  with  mobile  fuel  cells. 

The  objective  of  this  study  is  to  investigate  the  performance 
of  ATR  of  methane  over  a  sulfide  nickel  catalyst  on  a  gamma 
alumina  support  under  thermally  neutral  conditions  with  differ¬ 
ent  compositions  of  feedstock.  Modelling  work  is  also  carried 
out  to  study  further  the  conversion  behaviour  inside  the  reactor. 

2.  Experimental  study 

2.7.  Equipment 

The  schematic  layout  of  the  equipment  used  in  this  study  is 
shown  in  Fig.  1.  The  flowrates  of  air,  methane  and  water  are 
automatically  controlled  by  a  computer  according  to  set  values. 
The  fed  mixture  is  heated  by  either  an  electrical  heater  in  the 
start-up  stage  or  by  hot  refomate  gas  under  steady  operations  to 
ensure  that  water  is  in  vapour  phase  when  entering  the  reactor. 

The  reactor  is  essentially  a  stainless- steel  tube  (inner  diame¬ 
ter  =  30  mm,  length  =  200  mm)  that  is  filled  with  a  sulfide  nickel 
catalyst.  Reactor  temperature  is  measured  at  5  points  along  its 
centre  length  by  five  thermocouples.  The  outer  reactor  wall  and 
all  the  gas  pipes  are  covered  by  thermal  insulation  to  minimize 
heat  losses.  The  reactor  operates  under  atmospheric  pressure. 
The  reformate  sample  gas  is  channelled  to  a  system  of  analyz¬ 
ers,  at  which  its  compositions  is  measured. 

2.2.  Catalyst 

The  catalyst  was  a  commercial  sulfide  nickel  catalyst  Ni- 
0309S,  supported  on  gamma  alumina  and  was  supplied  by 
Engelhard.  This  is  a  new  type  of  commercial  catalyst  that 
is  specially  used  for  hydrogen  production  from  hydrocarbon 
fuel  reforming.  The  reasons  for  choosing  this  catalyst  is  that 
it  is  highly  active  and  relatively  cheap  and  has  prospective 
use  in  industries  compared  with  other  types  of  conventional 
nickel/alpha  alumina  catalyst.  The  catalyst  is  of  spherical  shape 
and  is  ready  for  use  as  supplied.  The  amount  of  catalyst  loaded 
in  the  reactor  was  68  g.  The  physical  properties  of  the  catalyst 
are  listed  in  Tables  la  and  b. 

2.3.  Experimental  procedure 

To  operate  the  reactor,  only  methane  at  a  flowrate  of  2 1  min- 1 
and  air  with  an  air-to-fuel  ratio  (A/F)  of  6  are  supplied  initially 
to  the  reactor  where  the  mixture  is  ignited  by  sparks  from  a 
spark-plug.  Heat  released  from  methane  combustion  heats  up 
the  catalyst  bed.  When  the  average  catalyst  temperature  reaches 
about  700  °C,  steam  is  added  to  the  inlet  mixture  and  then 
methane,  and  the  A/F  and  molar  water-to-fuel  (W/F)  ratios  are 
adjusted  to  the  desired  values  for  testing.  The  feedstock  condi¬ 
tions  are  varied  in  the  range  that  ensures  that  the  reactor  sustains 
a  thermally  neutral  operation  between  500  °C  and  900  °C.  Based 
on  a  theoretical  analysis  of  the  optimum  conditions  for  autother¬ 
mal  reforming  of  methane  [11],  the  A/F  is  varied  from  2.5  to  5 
and  the  W/F  is  varied  at  each  A/F.  The  W/F  is  increased  step- 
by-step  with  increments  of  0.5  from  a  starting  value  of  1  with 
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Fig.  1.  Schematic  layout  of  the  rig. 
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Table  la 


Catalyst  properties 


Nickel  content  (wt.%) 

9.8 

Sulfur  content  (wt.%) 

4.9 

Alumina  content  (wt.%) 

Balance 

Surface  area  (m2  g-1) 

155 

Total  pore  volume  (mg-1) 

0.9 

Size  of  the  sphere  (mm) 

1.75 

Average  crush  strength  (N) 

25 

a  time  interval  between  successive  steps  of  at  least  30  min  for 
the  reactor  to  reach  a  steady- state  condition.  The  reformate  gas 
concentration  is  recorded.  The  W/F  is  increased  until  the  reactor 
cannot  sustain  autothermal  reactions. 

2.4.  Experimental  results 

The  extensive  experimental  data  were  recorded  and  pro¬ 
cessed.  Typical  results  under  a  methane  flowrate  of  1201h_1 
are  expressed  in  Figs.  2-15.  These  show  the  effect  of  the  A/F 
the  W/F  on  the  performance  of  the  autothermal  reactor. 

The  variation  of  wet  and  dry  product  gas  concentrations 
against  W/F  at  A/F  =  3  are  presented  in  Figs.  2  and  3,  respec¬ 
tively.  Within  the  range  of  W/F  where  the  autothermal  reac¬ 
tions  is  sustained,  the  residual  methane  increases  (i.e.,  methane 
conversion  decreases)  a  little  with  increasing  W/F.  Carbon 


Table  lb 

Kinetic  parameters 


Reaction 

K0j  (molkgcat  1  s  !) 

Ej  ( J  mol  1 ) 

1 

9.287  x  102  bar-1'5 

30800 

2 

9.048  x  1011  bar0  5 

209500 

3 

5.43  x  105  bar-1 

70200 

4 

2.14  x  109  bar0  5 

211500 

monoxide  decreases  and  carbon  dioxide  increases  gradually  with 
increasing  W/F  due  to  the  improved  water-gas  shift  reaction. 
On  the  other  hand,  the  content  of  water  vapour  increases  dra¬ 
matically  in  the  wet  product  because  the  increase  in  the  steam 
reforming  reaction  is  not  directly  proportional  to  the  increase 
in  the  W/F  ratio.  As  a  result,  the  wet  hydrogen  concentration 
decreases,  although  its  concentration  increases  in  the  dry  prod¬ 
uct.  By  contrast,  the  nitrogen  concentration,  decreases  due  to 
increase  in  the  concentrations  of  water  vapour  and  carbon  diox¬ 
ide  in  the  product. 

The  methane  conversion  efficiency,  which  is  defined  as  the 
difference  in  the  moles  of  methane  between  the  inlet  and  out¬ 
let  of  the  reactor  divided  by  the  moles  of  methane  supplied, 
at  different  A/F  and  W/F.  The  results  show  that  methane  con¬ 
version  increases  with  increasing  A/F  from  2.5  to  3.5  and  then 
approaches  100%  in  the  full  range  of  W/F  under  study.  With 
A/F  <3.5,  methane  conversion  deceases  with  increasing  of  W/F. 


Fig.  2.  Wet  product  gas  concentrations  at  different  W/F  ratios  (A/F  =  3, 
VCh4  =  1201b-1). 


D.L.  Hoang  et  al.  /  Journal  of  Power  Sources  159  (2006)  1248-1257 


1251 


Fig.  3.  Dry  product  gas  concentrations  at  different  W/F  ratios  (A/F  =  3, 
yCH4  =  i20ih-1). 


Fig.  6.  Variation  of  hydrogen  wet  concentration  vs.  W/F  ratio  at  different  A/F 
ratios  (Vch4  =  1201h_1). 
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Fig.  4.  Methane  conversion  vs.  A/F  ratio  at  different  W/F  ratios  Fig.  7.  Variation  of  hydrogen  dry  concentration  vs.  W/F  ratio  at  different  A/F 
(VcH4  =  1201h  1  )•  ratios  (Vch4  =  1201  h_1). 


The  variation  of  hydrogen  concentration  in  wet  and  dry  prod¬ 
ucts  is  shown  in  Figs.  6  and  7,  respectively,  under  different 
A/F  and  W/F.  The  concentration  in  both  wet  and  dry  products 
decreases  with  increasing  of  A/F  due  to  the  increase  of  nitrogen 
in  the  products.  With  a  fixed  A/F,  the  wet  hydrogen  concen¬ 
tration  decreases  with  increasing  W/F  within  the  range  studied, 
whereby  the  reactor  can  sustain  autothermal  reactions.  This  does 
not  mean  that  the  rate  of  hydrogen  produced  decreases  with 


Fig.  5.  Methane  conversion  vs.  W/F  ratio  at  different  A/F  ratios 
(Vch4  =  1201h_1). 


increasing  W/F  because  of  the  increased  water  vapour  content 
in  the  products.  When  removing  the  water  vapour  from  the  prod¬ 
ucts,  the  dry  hydrogen  concentration  increases  with  increasing 
W/F,  as  indicated  in  Fig.  7. 

The  effect  of  W/F  on  the  flow  of  hydrogen  produced  can  be 
seen  more  clearly  in  Fig.  8,  which  shows  that  the  moles  of  hydro¬ 
gen  produced  per  mole  of  methane  fed  into  the  reactor  increases 


Fig.  8.  Mole  of  hydrogen  produced  per  mole  of  methane  fed  into  the  reactor  vs. 
W/F  ratio  at  different  A/F  ratios  (Vch4  =  1 20 1  h—  1 ). 
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Fig.  9.  Mole  of  hydrogen  produced  per  mole  of  methane  fed  into  the  reactor  vs. 
A/F  ratio  at  different  W/F  ratios  (Vch4  =  1 20 1  h“ 1 ). 


Fig.  10.  Variation  of  carbon  monoxide  wet  concentration  vs.  W/F  ratio  at  dif¬ 
ferent  A/F  ratios  (Vch4  =  1201h_1). 

with  increasing  W/F.  On  lowering  the  A/F,  the  maximum  W/F 
at  which  the  reactor  can  sustain  an  autothermal  reaction  is  low¬ 
ered  due  to  the  reduced  heat  release  to  vaporize  water  and  to 
compensate  the  endothermic  steam  reforming  reaction. 

The  effect  of  A/F  on  the  moles  of  H2  produced,  on  the  other 
hand,  is  different  as  can  be  seen  in  Fig.  9.  The  moles  of  H2  pro¬ 
duced  reaches  a  peak  for  A/F  ratios  between  3  and  3.5  with  a 
W/F  ratio  of  2-2.5.  Under  these  mixture  conditions,  the  maxi¬ 


Fig.  11.  Variation  of  carbon  monoxide  dry  concentration  vs.  W/F  ratio  at  dif¬ 
ferent  A/F  ratios  (Vch4  =  1201h_1). 


Fig.  12.  Mole  of  carbon  monoxide  produced  per  mole  of  methane  fed  into  the 
reactor  vs.  W/F  ratio  at  different  A/F  ratios  (Vch4  =  120  lh_1). 


Fig.  13.  Mole  of  carbon  monoxide  produced  per  mole  of  methane  fed  into  the 
reactor  vs.  A/F  ratio  at  different  W/F  ratios  (Vch4  =  1 20 1  h_1 ). 

mum  moles  of  hydrogen  is  about  1.8  moles  per  mole  of  methane 
supplied  and  the  wet  and  dry  hydrogen  concentrations  reach  28 
and  40%,  respectively. 

The  variation  of  CO  concentrations  and  CO  moles  against 
A/F  and  W/F  ratios  is  shown  in  Figs.  10-13.  Both  wet  and 
dry  CO  concentrations  and  the  moles  of  CO  produced  per  mole 
of  methane  dramatically  decrease  with  increasing  W/F.  This  is 


Fig.  14.  Catalyst  temperature  along  reactor  at  different  W/F  ratios  (A/F  =  3, 
VCh4  =  1201b-1). 
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Fig.  15.  Catalyst  temperature  along  reactor  at  different  A/F  ratios  (W/F  =  2, 
Vch4  =  120  lh1). 

because  of  increase  in  the  water-gas  shift  reaction  in  the  reactor. 
The  effect  of  A/F  on  CO  generation  is  similar  to  that  on  hydrogen 
generation.  With  an  A/F  of  3-3.5  and  a  W/F  of  2-2.5,  the  wet 
CO  concentration  is  as  low  as  5-7%,  the  dry  CO  concentration 
is  7-9%  and  the  moles  of  CO  is  as  low  as  0.4-0.45  per  mole  of 
methane  supplied. 

The  distribution  of  catalyst  temperature  in  the  reactor  at 
different  W/F  and  A/F  ratios  is  shown  in  Figs.  14  and  15,  respec¬ 
tively.  In  general,  the  catalyst  temperature  is  lowest  in  the  front 
part,  highest  in  the  middle  part  and  lower  in  the  rear  part  of 
the  reactor.  This  can  be  explained  as  follows:  the  front  part  is 
cooled  by  the  inlet  mixture  of  low  temperature,  the  middle  part 
is  heated  by  heat  release  from  the  total  combustion  of  methane, 
while  the  endothermic  steam  reforming  reaction  is  prominent, 
in  the  rear  part  the  reactor.  Increasing  W/F  leads  to  more  cool¬ 
ing  of  the  front  part,  which  causes  the  temperature  in  this  area 
to  decrease.  As  a  result,  more  of  the  total  combustion  reaction 
of  methane  shifts  to  the  middle  part  and  increases  temperature 
there  (Fig.  14).  At  a  fixed  W/F,  an  increase  in  A/F  means  an 
increase  in  the  combustion  rate  and  heat  release  and,  hence,  the 
reactor  temperature  increases  (Fig.  15). 

The  effect  of  inlet  methane  flowrate  on  the  reforming  per¬ 
formance  is  presented  in  Figs.  16-18,  which  show  gas  concen¬ 
trations  in  the  wet  product,  methane  conversion  and  moles  of 


Fig.  16.  Wet  product  gas  concentrations  at  different  methane  flowrates  (A/F  =  3, 
W/F  =  2). 


CH4  flowrate  (l/h) 

Fig.  17.  Methane  conversion  vs.  methane  flowrate  (A/F  =  3,  W/F  =  2). 

hydrogen  produced  per  mole  of  methane  under  different  flow 
rates  of  methane  feed  at  A/F  =  3  and  W/F  =  2.  Increasing  the 
rate  of  methane  feed  beyond  250 lh-1  will  decrease  methane 
conversion  (Fig.  17).  This  leads  to  a  decrease  in  H2  concentra¬ 
tion  in  the  products  and  the  moles  of  H2  produced  per  mole  of 
CH4  feed  (Fig.  18).  The  decrease  in  methane  conversion  with 
increasing  feed  flowrate  is  due  mainly  to  the  change  of  tem¬ 
perature  profile  in  the  reactor.  With  the  flowrate  of  feed  gas 
exceeding  a  certain  value  (250 1  h_1  in  the  present  experiment), 
a  higher  flowrate  will  cause  more  cooling  of  the  front  part  of  the 
reactor  and  result  in  more  inactive  catalyst  and,  hence,  a  lower 
chemical  reaction  rate. 

Based  on  the  results  presented  and  discussed  above,  it  can  be 
seen  that  the  optimum  operating  conditions  for  the  autother- 
mal  reactor  to  achieve  high  hydrogen  yield,  a  high  methane 
conversion  and  a  high  number  of  moles  of  hydrogen  produced 
per  mole  of  methane  feed  are  an  A/F  of  3-3.5  and  a  W/F  of 
2-2.5  under  a  flowrate  of  methane  of  120-250 1  h_1 .  With  these 
conditions,  methane  conversion  can  reach  95-99%;  to  the  hydro¬ 
gen  concentration  can  reach  39-41%  of  dry  products  (28-30% 
wet  products);  the  moles  of  hydrogen  can  reach  1.8  per  mole 
of  methane  feed  CO  concentration  can  be  as  low  as  7-9%  dry 
(5-7  %  wet)  and  the  number  of  CO  moles  can  be  as  low  as  0.4-0. 5 
per  mole  of  methane  feed.  Under  these  feed  conditions,  the  max¬ 
imum  catalyst  temperature  will  not  exceed  850  °C. 


Fig.  18.  Mole  of  hydrogen  produced  per  mole  of  methane  fed  into  the  reactor 
vs.  methane  flowrate  (A/F  =  3,  W/F  =  2). 
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3.  Mathematical  model 


3.1.  Chemical  reaction  scheme 


In  addition  to  the  experimental  investigation,  modelling  work 
has  been  conducted  to  examine  further  the  conversion  behaviour 
of  the  ATR  process  inside  the  reactor.  The  model  includes 
detailed  reactions  associated  with  partial  oxidation  and  steam 
reforming.  The  detailed  analysis  of  chemical  reactions  in  these 
processes  to  determine  the  reaction  scheme  for  ATR  has  been 
presented  elsewhere  [10,11].  The  main  reactions  considered  in 
the  model  are  as  follows: 


CH4  +  202 

->  C02  +  2H20  A//i(298)=  —802,  OOOkJkmol-1  (1) 

ch4  +  h2o 

CO  +  3H2  A//2(298)=  206,  000  kJ  kmol- 1  (2) 


CO  +  H20 

C02  +  H2  A//3(298)  =  — 41 ,  000  kJ kmol  *  (3) 


CH4  +  2H20 

C02  +  4H2  AH 4(29$)  =  165,  000  kJ  mol  *  (4) 


Thus,  the  model  takes  into  account  four  reactions  (1)— (4)  and 
six  gas  species,  i.e.,  methane,  oxygen,  carbon  dioxide,  water, 
carbon  monoxide  and  hydrogen,  nitrogen  present  in  the  inlet  air 
is  considered  as  a  diluent,  which  affects  only  the  gas  properties. 

Kinetic  rate  equations  and  kinetic  data  applied  for  sulfide 
nickel  catalyst  in  this  work  are  adopted  from  earlier  studies  with 
the  same  catalyst  [10].  The  rate  equations  and  kinetic  data  for 
Eqs.  (1)— (4)  are  as  follows: 


Ri  = 


klPCH4Pol 


(1  +  KCH4PCH4  +  K&pfc) 


C  1/2. 


(5) 


(pch4Ph2o  — 


I’L  Pco  \ 

Ke  2  J 


Table  2 


Equilibrium  constants 


Reaction 

Equilibrium  constant  Kej 

2 

5.75  x  1012  x  e-95411/(^r)  (bar2) 

3 

i.26  x  nr2  x  CS56S/(RT) 

4 

7.24  X  1010  X  q-^9960/{RT)  (bar2) 

Table  3 

Adsorption  constants 

Species 

Koi  (bar  l) 

A  Hi  (J  mol) 

CH4  (Combustion) 

2.02  X  10“3 

-36330 

02  (Combustion) 

7.4  x  10-5  bar0  5 

-57970 

ch4 

1.995  X  10“3 

-36650 

CO 

8.11  x  10“5 

-70230 

H2 

7.05  x  10“9 

-82550 

h2o 

1.68  x  104  bar 

85770 

reaction  j  (j  =  2-4)  Table  2;  Kf  =  x  q~ah?/rt  the  adsorp¬ 
tion  constant  of  species  i  (i  =  C H4,  02)  in  oxidation  reaction 
(1);  Ki  =  K0i  x  Q~AHi/RT  the  adsorption  constant  of  species  i 
(i  =  C O,  H2,  CH4,  H20)  in  reforming  reactions  (2)-(4),  which 
can  be  found  in  Table  3. 

The  rate  of  consumption  or  formation  of  an  individual  gas 
species  based  on  reactions  (1)— (4)  is  determined  by  summing 
up  the  reaction  rates  of  that  species  in  all  four  reactions,  i.e., 


O7H4  =  —R 1  -  R2  -  R4 

(9) 

ro2  =  -2R\ 

(10) 

rco2  =  R\  +  R3  +  ^4 

(11) 

m2o  =  2R\  -  R2  -  R3  -  2 R4 

(12) 

rco  =  R2  ~  R3 

(13) 

rn2  =  3  R2  +  R3  +  4  R4 

(14) 

where  r;  is  the  conversion  rate  of  gas  species  i  (i  =  CH4, 02,  etc.). 
3.2.  Model  development 


P  h2 


[pcoPh2o  ~ 


PH2PCQ2\ 
Ke  3  J 


li4  2 

~X5  Pch4Pu2o 
P  H2  V 


pf  PCO,  \  J_ 

Ke 4  J  X  Q? 


(7) 

(8) 


Qr  =  1  +  K CO  PCO  +  Kr2PH2  +  Kch4PCH4  + 


^h2o/?h2o 

^h2 


where  Rj  (kmolkgcat-1  h-1)  is  the  rate  of  reaction  j  (j=  1-4); 
pcn4>  Po2,  etc.  are,  respectively,  the  partial  pressures  of  gas 
species  CH4,  02,  etc.;  kj  =  k0j  x  e~EdRT  is  the  kinetic  rate 
constant  of  reactions  j  (j=  1-4)  and  is  determined  from  Ref. 
[10].  The  kinetic  data  are  shown  in  Tables  la  and  b,  where 
k0j  is  a  constant,  Ej  (kJkmol-1)  the  activation  energy;  R 
(kJkmol-1  K_1)  the  universal  gas  constant;  T  (K)  the  gas  tem¬ 
perature  in  the  reaction  zone;  Kej  the  equilibrium  constant  of 


The  schematic  layout  of  the  reactor  is  shown  in  Fig.  19.  The 
reactor  is  of  a  cylindrical  shape  (30  mm  in  diameter  and  200  mm 
in  length)  and  is  filled  with  a  spherical  sulfide  Ni  catalyst  with  a 
size  of  1.75  mm  and  9.8%  active  metal.  Radiation  in  a  catalyst 
bed  can  be  significant  only  at  temperatures  above  1000  °C  [17]. 
The  maximum  catalyst  temperature  in  this  study  is  less  than  this 


►  * 


Fig.  19.  Schematic  layout  of  2-D  reactor. 
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limit,  and  thus  radiative  heat  transfer  can  be  omitted.  Moreover, 
because  the  gas  flowrate  is  relatively  low  and  the  void  fraction 
of  the  bed  is  high  (0.35),  the  pressure  drop  is  assumed  to  be  neg¬ 
ligible.  To  simplify  the  model,  it  may  be  assumed  that  the  gas 
flow  in  the  reactor  is  uniform,  and  thus  a  1-D  model  can  be  used 
for  the  reactor.  The  radial  dispersion  of  the  gas  flow  due  to  a  cat¬ 
alyst  particle  size  of  1 .75  mm  and  the  radial  heat  transfer  due  to 
heat  loss  through  the  wall  of  the  reactor  may  have  some  effects 
on  the  overall  reforming  performance.  To  take  these  effects  into 
account,  a  2-D  reactor  model  is  more  appropriate  to  describe  the 
reforming  behaviour  under  the  mentioned  conditions  [11,18]. 
Hence,  a  2-D  model  is  developed  in  this  study.  Since  the  dis¬ 
persion  coefficient  in  the  axial  direction  is  much  smaller  than 
the  effect  of  gas  velocity,  it  can  be  ignored  without  significant 
influence  on  the  calculation  results.  In  addition,  the  heat  con¬ 
ductivity  of  the  gas  is  much  smaller  than  that  of  the  catalyst  bed, 
and  therefore  can  be  omitted.  In  the  gas  phase,  changes  in  the 
gas  properties  and  concentrations  with  time  are  normally  much 
smaller  than  those  with  space,  especially  when  the  operation  is 
in  a  quasi-steady  state;  hence  they  can  be  ignored  within  a  small 
step  time. 

With  the  assumptions  mentioned  above,  the  basic  governing 
equations  of  the  model  based  on  the  mass  and  energy  balance 
for  the  gas  phase  and  solid  phase  of  the  reactor  are  as  follows: 


where  dv  is  pellet  diameter;  Dt  the  gas  diffusivity  of  species  i  to 
the  mixture  of  the  other  gas  in  the  reactor;  Tbed  is  tortuosity  of 
the  bed  and  is  correlated  to  the  void  fraction  of  the  catalyst  bed 
s  as  follows  [20] : 

1 

Tied  — 

Based  on  the  phenomena  of  the  gas  flow  and  operating  con¬ 
ditions  of  the  reactor,  the  initial  and  boundary  conditions  are  set 
as  follows: 


initial  condition  t  =  0  :  T  =  T0; 

(18) 

at  the  reactor  inlet  face  z  =  0  :  Tg  = 

yin.  y 

yin 

— 

(19) 

dQ 

at  the  reactor  outlet  face  z  =  L  :  - 

dz 

=  0; 

II 

o 

(20) 

dCi 

at  the  reactor  centre  r  =  0  :  - =  0 

dr 

dTg 

>;  — ^ 

dr 

=  0 

(21) 

With  these  conditions,  Eqs.  (15a)  and  (17a)  can  be  replaced 
by  Eqs.  (15b)  and  (17b),  respectively,  i.e., 


T  PcatT 


(15b) 


dCi  f  d2Ci  1  dCi  \ 

M— —  —  Ddpi  (  2  - (T-  )  +  PoiO'i 

dz  \  drz  r  dr  J 


/ir  x  dT  (  d2T  d2T 

(15a)  Pb^pb tt-  —  K  (  + 


dt 


dr 2  dz 2 


dTg 

£Pgcpg 


dTg 

—  WPgCpg-^— 


+  Shh(T 


(16) 


4 

+  Sbh(Tg  —  T)  +  pcatJ2(-*Hj)Rj  (17b) 

7=1 


dT 

PbCpb—  =  K 
dt 


d2T  1  dT  d2T 
dr 2  r  dr  dz2 


4 

+  Shh(Tg  —  T)  +  pCSLt^^(—AHj)Rj 

7=1 


(17a) 


where  i  denotes  the  gas  species;  j  denotes  reaction  index;  pg,  pcat, 
Pb  (kg  m-3)  the  densities  of  gas,  catalyst  and  bulk  catalyst  bed, 
respectively;  cpg  and  cpb  (Jkg-1  K-1)  the  specific  heats  of  the 
gas  and  the  catalyst  bed,  respectively;  s  the  void  fraction  of  the 
catalyst  bed;  h^i  (ms-1)  the  mass  transfer  coefficient  of  the  gas 
component  i;  h  (W  m-3  K)  the  heat  transfer  coefficient;  T  and  7k 
(K)  the  temperature  of  solid  phase  and  gas  phase,  respectively; 
C[  (mol  m-3)  the  concentration  of  gas  species  i;  r  and  z  (m)  are 
cylindrical  coordinates;  5h  (m-2  m-3)  the  heat  transfer  area  per 
volume  of  the  catalyst  bed;  AHj  (J  mol-1)  the  heat  of  reaction 
j ;  K  (W  m-1  K-1)  the  heat  conduction  coefficient  of  the  catalyst 
bed;  the  dispersion  coefficient  of  gas  component  i;  u  (m  s- 1 ) 

is  the  superficial  gas  velocity  equal  to  the  ratio  of  the  volume 
flow  rate  to  the  cross-section  area  of  the  reactor. 

The  dispersion  coefficient  of  gas  in  a  catalyst  bed  is  dependent 
on  molecular  gas  diffusion,  bulk  gas  velocity  and  pellet  diameter 
and  can  be  expressed  as  [19]: 

(  Di 

Dfoi  =  s  I  - T-  0.5*7pw 

V  Tied 


At  the  interfacial  surface  of  the  inner  reactor  wall  and  the 
catalyst  bed  r  =  R : 


dCi 

dr 


=  0: 


(22) 


where  ra  is  the  ambient  temperature;  a  the  overall  heat  transfer 
coefficient  through  the  reactor  wall;  K  (W  m-1  K"  ])  is  the  heat 
conduction  coefficient  of  the  catalyst  bed. 

The  heat  transfer  coefficient  between  the  catalyst  and  the  gas, 
h  (W  m-2  K-1),  is  determined  using  the  Colburn  factor  [21]. 
The  overall  heat  transfer  coefficient  through  the  reactor  wall,  a , 
is  determined  from: 


1  _  1  b  1 

a  h{  X  hQ 


(23) 


where  b  is  thickness  of  the  reactor  wall  including  the  insula¬ 
tion  layer  (m);  hx  and  h0  are  the  heat  transfer  coefficient  on 
the  inside  and  outside  of  the  reactor  wall  (W  m-2  K-1),  respec¬ 
tively,  while  A  the  heat  conduction  coefficient  of  the  reactor  wall 
(Wm-1  K-1),  which  can  be  taken  from  Ref.  [22];  h[  is  given 
by  Cussler  [23]  and  hQ  for  free  convection  is  taken  from  Ref. 
[11,24]. 

The  set  of  three  governing  Eqs.  (15)— (17)  with  Eqs.  (18)-(22) 
for  the  initial  and  boundary  conditions,  combined  with  the  heat 
and  mass  transfer  coefficients  is  then  solved  for  the  temperature 
and  gas  concentration  along  and  across  the  reactor  using  a  finite 
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Fig.  20.  Comparison  between  predicted  and  measured  methane  conversions  at 
different  W/F  and  A/F  ratios  (Vch4  =  1201  h-1 ). 


Relative  length  of  the  reactor 


Fig.  22.  Predicted  catalyst  temperature  along  the  reactor  at  different  W/F  ratios 
(A/F  =  3,  Vch4  =  1201h_1). 


difference  method.  The  overall  temperature  and  compositions  of 
the  products  are  obtained  from  the  predicted  values  at  different 
space  elements  at  the  rear  face  of  the  reactor. 

3.3.  Modelling  results 

The  program  code  is  run  under  initial  catalyst  temperature 
Tq  =  673  K  (400  °C),  at  which  the  reactor  can  light-off  and  sus¬ 
tain  auto  thermal  reactions  [11].  The  inlet  mixture  compositions 
and  flowrate  are  set  to  the  same  values  as  that  of  the  experiment 
carried  out  previously.  When  the  reactor  is  at  a  fully  steady- state 
(all  output  parameters  are  stable),  the  gas  parameter  inside  the 
reactor  and  the  product  compositions  are  extracted.  Some  typi¬ 
cal  predicted  data  under  a  fuel  flowrate  of  Vch4  = 120 lh-1  are 
presented  in  Figs.  20-25. 

The  predicted  and  measured  conversion  of  methane  under  dif¬ 
ferent  W/F  and  A/F  are  compaired  in  Fig.  20.  There  is  quite  good 
agreement  between  the  two  sets  of  data.  The  agreement  between 
modelling  and  experimental  results  is  also  seen  in  Fig.  21,  which 
shows  the  dry  gas  product  concentrations  at  A/F  =  3  and  differ¬ 
ent  W/F.  With  this  successful  validation,  the  simulation  program 
code  can  be  used  to  study  the  conversion  behaviour  inside  the 
reactor. 


Relative  length  of  the  reactor 

Fig.  23.  Predicted  methane  conversion  inside  the  reactor  at  different  W/F  ratios 
(A/F  =  3,  Vch4  =  1201h_1). 

The  catalyst  temperature  along  the  reactor  at  A/F  =  3  and  dif¬ 
ferent  W/F  is  given  in  Fig.  2.  Increasing  water  content  at  the  inlet 
leads  to  a  change  in  the  temperature  profile  of  the  catalyst  bed. 
The  front  part  of  the  reactor  is  cooled  down  with  increasing  W/F. 
This  then  causes  the  methane  conversion  to  decrease  sharply  in 
the  front  part  of  the  reactor,  as  seen  in  Fig.  23,  and  hence  results 
in  a  decrease  in  the  overall  conversion. 

The  formation  of  hydrogen  and  carbon  monoxide  in  the  reac¬ 
tor  is  shown  in  Figs.  24  and  25,  respectively,  the  decrease  in 


Fig.  21.  Comparison  between  predicted  and  measured  gas  concentrations  at 
different  W/F  ratios  (A/F  =  3,  Vch4  =  120 1  h_1). 


Fig.  24.  Predicted  wet  hydrogen  concentration  inside  the  reactor  at  different 
W/F  ratios  (A/F  =  3,  VCH4  =  ^Olh-1). 
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Relative  length  of  the  reactor 

Fig.  25.  Predicted  wet  carbon  monoxide  concentration  inside  the  reactor  at  dif¬ 
ferent  W/F  ratios  (A/F  =  3,  Vch4  =  1201h_1). 

hydrogen  formation  with  increasing  W/F  is  a  consequence  of 
the  decrease  in  methane  conversion  presented  in  Fig.  23.  The 
decrease  of  carbon  monoxide,  on  the  other  hand,  is  due  to  the 
improved  water-gas  shift  reaction  at  lower  temperature  and  high 
water  content. 

The  simulation  code  allows  determination  of  all  the  gas 
parameters  and  the  catalyst  temperature  at  any  position  inside 
the  reactor  and,  hence,  it  can  be  used  to  study  the  factors  that 
influence  the  reforming  performance  and  for  optimum  reactor 
design. 

4.  Conclusions 

Autothermal  reforming  of  methane  over  a  sulfide  nickel  cata¬ 
lyst  is  investigated  experimentally  and  theoretically.  The  exper¬ 
iments  are  carried  out  under  thermally  neutral  conditions  in  a 
cylindrical  reactor  of  30  mm  in  diameter  and  200  mm  in  length. 
The  results  show  that  the  conversion  behaviour  of  the  reactor 
strongly  depends  on  A/F,  W/F  and  the  inlet  mixture  flowrate. 
The  optimum  condition  of  autothermal  methane  reforming  is 
found  at  a  molar  A/F  of  3-3.5,  a  molar  W/F  of  2-2.5  and  a 
fuel  flowrate  below  250 1  h_1 .  Under  these  conditions,  methane 
conversion  is  95-99%;  the  hydrogen  yield  the  39-41%  on  a  dry 
basis,  and  1  mole  of  methane  can  produce  1.8  moles  of  hydro¬ 
gen.  The  carbon  monoxide  concentration  is  as  low  as  7-9%  on 
a  dry  basis  (5-7%  on  a  wet  basis),  or  0.4-0. 5  mole  per  mole  of 
methane  feed. 


A  two-dimensional  heterogeneous  reactor  model  has  been 
successfully  developed  and  validated.  The  study  considers  four 
main  simultaneous  chemical  kinetic  reactions  that  involve  six 
species  with  detailed  kinetic  conversion  and  heat  and  mass  trans¬ 
fer  phenomena  in  the  reactor.  The  model  is  useful  for  analysis 
of  the  factors  that  influence  the  reforming  performance  and  for 
optimum  reactor  design. 
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